Symplectic parallel integrators in the realm of 
Hofer's geometry 

Hugo Jimenez-Perez 
January 4, 2012 



Abstract 

Symplectic integrators constructed from Hamiltonian and Lie formalisms 
are obtained as symplectic (indeed Hamiltonian) maps whose flow follows 
the exact solution of a "sourrounded" Hamiltonian H — H + h k Hi. Those 
modified Hamiltonians depends virtually on the time by the timestep size 
h. When the numerical integration of a Hamiltonian system involves 
more than one symplectic scheme as in the parallel-in-time algorithms 
and specifically the Parareal scheme, there are not a simple way to con- 
trol the dynamical behavior of the error Hamiltonian. The interplay of 
to different symplectic integrators can degenerate their behavior if both 
have different dynamical properties, reflected in the number of iterations 
to have a good approximation to the final sequential solution. Considered 
as flows of time-dependent Hamiltonians we use the Hofer's geometry to 
search for the optimal coupling of symplectic schemes. As a result, we 
obtain the constraints in the Parareal method to have a good behavior for 
Hamiltonian dynamics. 

1 Introduction 

Symplectic integrators are the natural methods for simulating Hamiltonian dy- 
namics. To construct a symplectic integrator, we can follow two different pro- 
cedures: on one side we can use the Hamiltonian formalism using generat- 
ing functions, Lie transforms, etc. On the other side we can take a general 
method, for instance, Runge-Kutta formulas, modifying the coefficients to sat- 
isfy the symplecticity conditions 1 . In general, we obtain implicit methods 
but when the Hamiltonian can be separated into kinetic and potential energies 
H{q,p) = T{p) — V(q) we can construct explicit methods easy to implement. 

They are constructed using the diffeomorphisms which left invariant the sym- 
plectic form u = dp A dq which defines de Hamiltonian vector field Xh ■ Such a 
diffeomorphisms are called symplectic diffeomorphisms or symplectomorphisms 
and they form a subgroup denoted by Symp(Af , u>). In particular, the flow of any 
Hamiltonian vector field is a symplectic diffcomorphism and all of them form 
another subgroup of diffeomorphisms called the Hamiltonian diffeomorphisms, 
denoted by Ham(M, ui). 
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Moreover, using the Lie formalism we can consider Symp(Afw) and Ham(M, ui) 
as Lie groups and the sets of symplectic and Hamiltonian vector fields as their 
Lie algebras at the identity element. Then, for every Hamiltonian vector field 
Xh £ f)am(M, w), its flow is given by the exponential (ft{xo) = e tXll x . For 
a fixed t — h G R + the mapping x(t ) h4 x(t + h) generated by iph = e hx " 
is a Hamiltonian map which defines the symplectic integrator. In fact, this is 
a Hamiltonian integrator which preserves more structure than the symplectic 
one. 

Classically the paralelization of this type of systems is performed by a decom- 
position of the phase space (domain decomposition) or looking for parallelizable 
tasks into the method or into the equations. However, in the last two decades 
there were several attempts to develop another type of parallelization for sym- 
plectic integrators based on the decomposition of the time variable. The first 
attempt of some parallel-in-time algorithm for a scalar differential equation was 
published by Nievergelt in 1964 [35]. The idea is to decompose the total time 
in several subintervals which can be computed in parallel. Each interval, called 
a "branch" , must be modified propagating the corrected local initial condition 
to each subinterval; this technique has derived in the multishooting methods. 
Although the Nievergelt 's algorithm is not iterative, almost all other algorithms 
use an iterative process to approximate the sequential numerical solution. Those 
iterative algorithms consist in two steps: one predictor which estimates in par- 
allel the value of several branches and one correcto$\ which approaches the 
final solution propagating the predictions between different branches. We call 
them the time-parallel algorithms and all of them differ in the corrector step 
which uses different iterative process to convergence. Of course, there are others 
differences but in this paper we are interested in the corrector step. 

As noted by Saha, Stadel and Tremaine [22] one way to time-parallelize an 
almost integrable Hamiltonian system is to compute in parallel several branches 
saving the perturbing contributions and to propagate them computing in se- 
quence the integrable Hamiltonian part. This technique is reproduced in [25] 
for high-order symplectic integrators however, it works fine if the ratio of the 
computing time of the integrable over the perturbing part is very small. This 
approach is very accurate but expensive. On the other extreme, there is the 
parareal method introduced by Lions, Maday and Turinici in |33j and refined in 
[4] . In this method the propagation of the predictions is made by a simple incre- 
ment of the corrector at every iteration. As a result, it is a very fast algorithm, 
however, for Hamiltonian systems there are several inconvenients associated to 
the non preservation of the geometric structure of the underlying integrators. 

In order to deal with this type of problems, Bal and Wu [5| have done the first 
step considering a new way to spread the information between branches in the 
sequential step and practically destroying the "pure" parareal scheme. Also, 
Dai et al. [T2] have introduced another variation of the parareal step, using 
symmetries and projections into the energy manifold to preserve the geometric 

1 Many authors inverse the terminology considering the parallel step as the corrector and 
the sequential step as the predictor. Our choise is evident when we relate the corrector step 
with the symplectic correctors studied in 50 , 37 38 
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properties of the underlying integrators. Recently, the author has proposed a 
geometric corrector step using Lie's algebras [28] which is equivalent to that 
from Dai et. al.. In this paper we translate that approach to the Hofer's 
geometry in order to search not only for Hamiltonian maps but for optimizing 
the energy (which we relate with the number of iterations) to go from the first 
guess solution to the final solution. 

2 Lie algebras and Hamiltonian vector fields 

We consider the phase space of a Hamiltonian system as a symplectic manifold^ 
M = T*M" = M 2 " with the canonical symplectic form w = dp A dq. Denote by 
X(M) the set of all vector fields and by $(M) = C°°(M) the set of all smooth 
functions over M. 

We define de binary operation [-, ■] : X(M) x X(M) — > X(M) by the rule 

[X,Y] = XY-YX, X,Y£X(M), (1) 

called the Lie bracket, which is: bilinear, alternating, and satisfies the Jacobi 
identity. The set of all the vector fields X(M) equipped with the Lie bracket 
(J!) obtains the structure of Lie algebra (X(M), [•,•]). 

Let CxF be the Lie derivative of F along the vector field X £ X{M). CxF 
meassures the change of F along X where F can be a function, a vector field, a 
p-form or, in general, a tensor. 

The Lie derivative of the symplectic form uj along X is given by the Cartan's 
magic formula 

Cxu — dixu + ixdu, (2) 

where d is the exterior differential and ix^J = oj(X, •) is the contraction of u> by 
X or equivalently the inner product of the vector field X with the 2-form uj. 

We say that the vector field X is symplectic if its flow preserves uj, which 
means £iw = 0. Since ui is the canonical symplectic 2-form then uj = —d\ 
where A = pdq is the Liouville form. Consequently, the second term in the right 
hand side of © is zero. In other words, a vector field X is symplectic if the 
1-form ix<-J is closed which means that it belongs to the kernel of d 

d(i X Lu) =d(w{X,-)) = 0. 

We denote the set of symplectic vector fields on M by sp(M, uj). 

We say that X is Hamiltonian if, in addition, ixoj is exact, i. e., there exists 
/ £ 5(M 2n ) such that 

l X LU = L0(X,-) = ~df. (3) 

We call / a Hamiltonian function for X and we write X = Xf to specify that X 
is the Hamiltonian vector field associated to /. the set of all Hamiltonian vector 

2 A11 the computations and results listed here apply to arbitrary symplectic manifolds 
(M,w). 
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fields on M is denoted by t)atn(M,uj). Finally, the triplet (M,uj,Xh) defines a 
Hamiltonian system over M. For general mechanical systems the configuration 
space can be consider as a Riemannian manifold (N, g) and the phase space 
becomes the cotangent bundle M = T*N which has a natural structure of 
symplectic manifold. 

We need a non obvious result from the theory of differential p-forms which 
give us the expression of the inner product of the Lie bracket with a p-form a 
over any differential manifold M. The inner product i\ x ,Y] a f° r anv p-form a 
is given by ((7J pp 73]): 

i[x,Y] a = Cxiiyoi) - iy(jO-xa). 
Consequently, for every two X, Y G sp(M, cu) we have 

i[x,Y]U = Cx(iYu)+Q 

— d (i x (iyw)) + 

= d{u(Y,X)) 

= -d(u{X,Y)) 

where we used £yw = and du) = 0. 

These computations has important consequences: 1) [X, Y] is a Hamiltonian 
vector field. 2) Since t)<xm(M, uS) C sp(M,u>), the Lie bracket [■, ■] gives them 
the structure of Lie subalgebras of X(M), and 3) t)am(Af, uS) is an (in fact the 
maximal) ideal of Sp(M, w) with respect to [•, •], i.e. 

[sp{M,u),sp(M,u)\ C \)am(M,w). 

We have the relations 

f)om(M, u) c sp(M, uj) c X(M). (4) 

It is a well-known fact that \)am(M. uj) = sp(M, to) if and only if the fundamental 
group of M is trivial; in other words, when M is simply connected. 

Now we link this point of view with the classical development in local coordi- 
nates. Let H : M — > R be a differentiable function with Hamiltonian vector field 
Xh on M. Select a point v G 4>^ 1 (M) in a local chart of M. The Darboux's 
theorem says that, locally, all symplectic manifolds are symplectomorphic to 
T*R ra = R 2 ™ and then we can consider that v G M 2 ™. In canonical symplectic 
coordinates v = (q,p) such that q G R" and p G T*R n = M™. The vector field 
Xh in local coordinates is 

. OH . dH 

which are called the Hamilton equations. 

For two differentiable functions f,g G S'(M), their associated Hamiltonian 
vector fields Xf,X g fullfils u)(Xf,X g ) = —{f,g} where {■,■} is the Poisson 
bracket for functions defined by 

{/ } ^ ^ 5 ^ 5 (6) 

dq dp dq dp 
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The binary operation © is bilinear, antisymmetric and fulfills the Jacobi iden- 
tity. The space of real- valued differentiable functions $(M), equipped with the 
Poisson bracket ©, obtains the structure of Lie algebra. It is possible to write 
the Hamiltonian vector field in terms of the Poisson brackets by 

z = X H (z) = {z,H}. 

There exists a natural anti-morphism of Lie algebras between the algebra of 
differentiable functions 3(M, {•, •}) and the algebra of Hamiltonian vector fields 
X(M, [;■}) given by 

/ h- X f (7) 
{/,</} * X {lg} = -[X f ,X g ] (8) 

Since the Poisson bracket of two functions is a function then the Lie bracket of 
two Hamiltonian vector fields is again a Hamiltonian vector field, as we have 
shown before. 



3 Hamiltonian diffeomorphisms and Hofer's ge- 
ometry 

A symplectic diffeomorphism or symplectomorphism of a symplectic manifold 
(M,ui) is a C°° diffeomorphism <j> £ Diff(M) which preserves the symplectic 
structure lo, it means that the pull-back of (f> fulfills 4>*(lo) = lo. The support 
of a diffeomorphism <fr is the closure of {x £ M\4>{x) ^ x}. The set of all 
symplcctomorphisms with compact support form a group, denoted Symp(M, u)) 
(with the law of composition of mappings). 

Let Xh € f)am(Af, lo) be a Hamiltonian vector field. The flow ipf — 
e tXH (xo) of Xh is a one parameter subgroup of symplectic diffeomorphismf[f|. 
The set of all the symplectomorphisms which arise as the flow ip{ of Hamil- 
tonian vector fields Xf form another subgroup of diffeomorphisms called the 
Hamiltonian diffeomorphisms, denoted by Ham (A/, lo). It is easy to prove that 
Ham(Af, lo) is an infinite-dimensional subgroup since for every / £ $(M), we 
have a mapping X_ : $(M) — > X(M) which maps / i-> Xf with kei(X_) = M 
the constant functions. The exponential map e - : f)om(Af, lo) — > Ham(Af, lo) is 
injective then the composition e - o X_(H) = e tx " sends an infinite dimensional 
basis of $(M) to an infinite-dimensional basis of Ham(M, lo). 

We can relate the groups Ham, Symp and Diff with the Lie algebras in Q 
by the exponential map as follows 

Ham(A4» C Symp(M,w) c Diff(M) 

exp | exp t exp t (9) 

t)am(M,co) c sp(M,u) c X(M). 

3 We use the exponential map of vector fields since Symp(M, w) is actually a Lie group 
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The reader must note that the exponential mapping is, in general, not surjec- 
tive and the group Symp(M, u>) can have several components. By construction 
Ham(M, u) belongs to the identity component of Symp(M, to). 

Consider the set of time-dependent Hamiltonian functions H : M x I — )• M 
with compact support. We can normalize such a functions since for every in- 
terval / = [0,a],a < oo the flow of the vector field associated to the function 
H(x, t) — aH(x, at) defined on / = [0, 1] is again a Hamiltonian flow. More gen- 
erally for every smooth function g(t) with <?(0) = the flow <Pg(t\ 1S Hamiltonian 

with normalized Hamiltonian function ^(t)H(x, g(t)) |41) . These properties of 
rescaling in time are used to regularize singularities in mechanical systems. 

We can restate the definition of Hamiltonian diffcomorphism in the following 
way: a Hamiltonian diffeomorphism is a diffeomorphism (j) : M — » M which can 
be written as the time-l-map of a Hamiltonian flow, i.e., <fi — ip^ for some 
time-periodic Hamiltonian H : S 1 x M — > R [45]. Let us denote by Ht the 
function H (t, •) on M . In the following, we will normalize the Hamiltonians and 
consider only the time-l-maps in Ham(M, w). We denote by T-L(M, to) the group 
of time-l-map of Hamiltonian flows which is a subgroup of Ham(M, w). 

Remark 1. There is a constraint in all these definitions since the theory ap- 
plies for functions and diffeomorphisms with compact support. However, for 
symplectic integrators, we do not need global properties and the most important 
thing is the numerical tests for the error behavior. 

A diffeomorphism <\> E Diff(M) is said to be isotopic to the identity if there 
exists a smooth map H : M x [0, 1] — ► M such that if ht : M —> M is given by 
h t (x) = H(x,t), then ht is a C°° diffeomorphism, ho — idu and h\ — </>. We 
say that h t is an isotopy from </> to the identity. 

We say that ht is a Hamiltonian isotopy if there exists a smooth family of 
functions H t : M — > ffi. such that 



In this way, we have constructed curves or trayectories in Ham(Af , lu) which 
connects any Hamiltonian diffeomorfism in T-L{M, u) with the identity map. 

For every <\> G Ham(Af, uS), choose a Hamiltonian isotopy $ = ((f> t ) from <f> to 
the identity. Hofer [51] defined the length of this isotopy by 



where osc(-fft) := max(_ff f ) — min(iJ t ) denotes the oscillation of a function on 
M. For H S H(M,u), it is clear that 1(H) = if, and only if, H = 0. The 
distance from the identity, or energy, of an element <f> G Ham(Af , lu) is defined 
as 



ix ht u = -dH t 



(10) 




(11) 



d(id,cj)) :=M{l(H)\H <EH,cj) = >p? }■ 



(12) 



G 



Let us extend the distance to a function d : Ham(M, uS) x Ham(M, uS) — > [0, oo) 
by setting d((j), tp) := d(id, tp o Hofer has showed in [23] that •) is a bi- 

invariant metric on Ham(Af, ui) defined intrinsecally. Then the set Ham(M, uj) 
with the metric (|12l) is called the Hofer's geometry and it is a fundamental stone 
in symplectic topology. 

As was pointed out by Siburg in [45 , to every Hamiltonian dynamical system 
corresponds one single path in Ham(M, uj) and vice versa. All the dynamical 
properties of the Hamiltonian system, including its periodic orbits, heteroclinic 
connections, etc. are contained in the Hamiltonian isotopy. Moreover, Baily 
and Poltcrivich have showed in [5] that the bifurcation diagram of every Hamil- 
tonian system is preserved for every Hamiltonian diffcomorphism contained in 
a geodesic path in Ham(M, uj). 

For instance, let H t 6 H(M,uj) be an admissible Hamiltonian which means 
that H t has a compact support for every t £ I. Then H t is said to generate a 
minimal geodesic if d(id 7 cpf) — 1(H). However, it is very difficult to work with 
the space of geodesies in Ham(Af, w). 

4 Symplectic integrators and Hamiltonian maps 

Symplectic integrators are the natural numerical methods for simulating Hamil- 
tonian dynamics. From the geometrical point of view, the most natural inte- 
grators are obtained by the Hamiltonian formalism and the Lie theory applied 
to the group Ham(M, uj) and its Lie algebra fjam(M, uj) = Tjd (Ham(A/,w)). 

Consider the Hamiltonian system (M,uj,Xh). The flow generated by the 
Hamiltonian vector field Xh £ t)Qm(M, uj) is the one-parameter subgroup of 
Ham(M, w) defined by <pf(xo) = e tXn (x ) where xo £ M is the initial condition 
of the vector field. For a fixed value t = t the mapping e TXn : M — > M is a 
symplectic, actually a Hamiltonian map. In the generic case, e TX " is very 
complicated and is given in (analytical) implicit form. 

However, there exists an important class of Hamiltonian functions for which 
the Legendre condition^ is satisfied. This class is formed by separable functions 
H{q,p) = T[p) + V(q) where T(p) is the kinetic energy and V(q) is the poten- 
tial. (We use the plus sign in order to develop the exponential as a product.) 
Then the Hamiltonian vector field is separable and each part can be integrated 
independently obtaining an implicit method. If we write Xh — Xt + Xy the 
flow generated by Xh becomes 

<pf = e TX " = e *-(*T+*v) 

and since each part can be integrated independently, we can estimate e rXT and 
e rXv directly. However, in general [Xt, AV] ^ and therefore 

e r(X T +X v ) ^ e rX Te rX v 

4 The Legendre condition ask for the convexity of the Hamiltonian function for which its 
Hessian does not vanish. 
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concide only in the first term. This becomes a first order method. 

In order to obtain higher order methods, we search for coefficients {0^}™ and 
{frill™ sucn that J2i a i = 1 an d J2i k == 1 to estimate e rJfH by the composition 
of maps 

e TX T +rX v ^ ^ciitXt ^>\tX v e a m rX T e a m rX v (^3) 

We must impose additional conditions to the coefficients a.; and hi for matching 
more terms on both sides of the expression (TT3|) . For that, we use the Baker- 
Campbell-Hausdorff (BCH) formula and derive a set of polynomial conditions 
for ai and hi. In particular, the derivation of time-parallel methods impose a 
reversibility condition which requires that the symplectic method be symmetric 

Using the BCH formula we can write the symplectic integrator (fT3|) of order 
k with its residue as 



e 



tXt+tXv _ gOi-rJCT gbirXv e a m rX T e a m rX v _|_ Q( T k + 1 ) ^4) 



which implies 

^tXt^tXv _ _ _ e o m rX re o m TX v = g rX„ _ 0(^+1), ( 15 ) 



The numerical solution follows a modified or "surrounded" Hamiltonian which 
we consider as the nonautonomous function 

H(t, q,p) = H(q,p) + t^H^t, q,p). (16) 

This property give us a measure for the error commitcd by the symplectic meth- 
ods which corresponds to the error Hamiltonian r fc+1 H\{t, q,p) for fixed r. The 
interested readers will find a deeper discussion on the subject in the works 
of Suzuki [47], Yoshida [51] for this type of Hamiltonians and McLachlan [37] . 
Laskar and Robutel [31] for Hamiltonians of the form H(q,p) = T(p) + eV(q) for 
small e. A more general discussion on the construction of symplectic integrators 
using the BCH formula is found in Hairer et.al. [22]. 

For every symplectic integrator obtained with this process, we have a Hamil- 
tonian function (fl6|) characterized by its error Hamiltonian which depends on 
({<!«}, {&i},r), and then an analytical flow (p H G Ham(M, u>) associated to the 
numerical solution. For f,g £ $o (M) such that / ^ g the Hamiltonian maps 

tpf =e rX f and <p g = e TX °, <p f , (p 9 € H(M , uj) , (17) 

associated to their flows fulfills <p* ^ f 9 ■ The Hofer's metric give us a way 
to measure the distance between them in Ham(M, uj). This fact permit us to 
compare symplectic integrators in an geometrical framework. 



5 Time-parallel integrators 

The corner stone in the theory of time-parallel integrators is the use of two 
different numerical flows, generally using a two level discretization 8t -< At -< 
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[0, T] and a corrector step which propagates sequentially with the coarse flow, 
the values obtained in parallel with the fine flow. In the seminal article of 
Nievergelt he proposed, as an example, using the Euler method on each 
discretization. This produces two numerical flows which interplays to estimate 
the final sequential solution. 

The first documented proposal for a time-parallel algorithm for Hamiltonian 
systems (using a symplectic integrator), is in the work of Saha, Stadel and 
Tremain [H]. They propose a time-parallel algorithm for almost integrable 
Hamiltonian systems in action-angle coordinates. In such type of coordinates, 
the integrable part of the system corresponds to the actions which are constants 
in the flow. Instead of having a two level discretization, they use a single level 
discretization with two flows: one for the integrable system and the other one 
for the complete (perturbed) system using the symplectic mid-point rule. In 
this work, at least from the examples they show, it is not possible to use a high 
order symplectic method. 

The next interesting proposal is the Parareal algorithm which was used for 
Hamiltonian systems with a lot of degrees of freedom [3J E]. However, for 
long time simulations of systems with a few degrees of freedom Parareal has not 
a well behavior [T31 [IS]. The problem is that the parareal algorithm (in fact, 
the parareal step) does not preserves the symplecticity when it propagates the 
fine flow with the coarse one. There are several attempts to obtain a better 
corrector step in order to preserve the symplecticity [SJ [TH [25] , however, there 
are not a concrete answer to this problem. Here we ask for an additional point. 

Given two different symplectic maps, one for each discretization we can 
obtain the symplectic map of its Lie bracket and construct a symplectic corrector 
as in [5, 28 . At each iteration corresponds a point in Ham(M, uS) and then there 
exists a Hamiltonian isotopy h t C Ham(M, w) which contains all such points. 
Since the fine flow is given (it corresponds to the final numerical solution we 
expect to approximate), the question rests in the choice of the coarse flow in 
order to approach the fine solution in the minimal number of iterations. We 
give a partial answer in the next paragraphs. 

Let's consider the initial value problem 

y(t) = X H (y(t)), y(0) = 2/o. (18) 

where y : [0, T] —> M and Xh € (jam(M, ui). We discretize the problem by 
partitioning [0, T] in N subintervals of size At = T/N which we call branches 
as in We set to — 0, t n — T and ti = iAt such that 

= t < h < t 2 < ■ ■ ■ < U < ■ ■ ■ < t N = T. (19) 

We write At -< [0, T] for this partition. Each branch in At is decomposed in 
Ns G N subintervals St = At/Ng which corresponds to the final resolution. Then 
St -< At -< [0,T]. 

On this discretization we introduce two levels of resolution: a coarse sym- 
plectic solver G on the At resolution and a fine symplectic solver F on the St 
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one. Their flows ip G and (p F arc defined uniquely by their coefficients and their 
timesteps ({di}, {bi}, At) and ({cj}, {di}, 5t) respectivelly. 

The fine and the coarse propagation of the solutions y n on a branch are 
respectively given by F(y n ) and G(y n )- For the problem (TT5)) we estimate the 

The 

(20) 
(21) 



first guess sequence with the coarse solver = y(0) and y^+\ = G{yn')- 
time-parallel algorithms are given by some variation of the to steps 



y ( o k+1) 



- y ( k \ 



where the subscripts are the propagation in time and the superscripts are the 
iterations. In (|21[) . T is the corrector which is an operator on the symplectic 
solvers. The general algorithm is as follows 



Time-parallel algorithm 1. 



1 

2 
3 
4 

5 

G 
7 

8 

9 
10 



Setup of the initial guess sequence 

y% = y(o),y° n+ i = G{yl) 

For fc = 1 to k max 

Parallel resolution on [T n ,T n+1 ]: 

compute T > 
For n = k to N 

Sequential corrections: 

compute yli+i^ = r (^F 
end for (n). 
end for (k). 



(fc) 

Vn 



(AH 



Vn' 



In particular the parareal algorithm implements (|21l) as 



y n +i 



Q[y, 



Q[yi k) 



(22) 



Expression (I22[) has been called the parareal iteration. 

Recall that the construction of a symplectic corrector uses the Lie bracket 
of Xp and Xq with at least one reversal integration, then we impose that both 
symplectic solvers be symmetrical |28j . 

Every implicit symplectic integrator for separable Hamiltonian systems can 
be uniquely determined by its coefficients {c{\ and {di}. Its numerical flow ip H 
by the coefficients and the timestep ({c,}, {di}, r). Then, we identify the flow of 
every symplectic integrator with timestep r with the triplet <p H = ({cj}, {di}, r). 
Suppose that we have, as the fine symplectic scheme of order fc, the more accu- 
rated integrator in the family of the fc-order symplectic integrators. 

Theorem 1. Given the two level discretization St -< At -< [0,T] and the final 
(optimal) symplectic integrator ({cj}, {di}, St) with flow ip F the closest flow ip G 
in Ham(M, oj) for At is given by ({cj}, {di}, At). 
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Proof. Immediate using the triangle's inequality property of the Hofer's 
metric. □ 

As a consequence, the use of some symplified symplectic scheme increases 
the number of iterations and, in general, introduces an erratic behavior since 
the dynamics of both flows is in general not equivalent. 

In this case, the corrector step constructed by the author in 28 coincides 
with those studied by Wisdom and Holman [50] McLachlan 36] and Laskar and 
Robutel [31]. 

There are more results from the symplectic topology which apply for partic- 
ular families of Hamiltonians. For example, a function H G $o(M,u)) is called 
to have quadratic growth if there exists c £ M finite such that \d 2 H{z)\ < c for 
all z € M. For such functions, the Hamiltonian isotopy of the Theorem [T] is a 
geodesic in Ham(Af,w) [35 , Ch. 12]. 

6 A numerical test 

We compare the pure parareal with the more accurate JL11 algorithm intro- 
duced in |26j for the one dimensional Spin-orbit problem with several values of 
the parameter a 

H(q,p) = ^p 2 -e(cos{2q) + a(cos2q + 9-7cos2q-9)) (23) 

We select a 8th order symplectic integrator from the SBAB n family for both T 
and Q solvers. 

The first approximation with the coarse solver approach faster than the 
symplectic proposed in [25J and JL11 proposed in [33]. However, as was showed 
in the later document, JL11 obtain exactly the same sequential solution. 
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Figure 1: Numerical test of the Spin-orbit problem with JL11 and parareal. 
The parammeters are e = 0.1, a = 0.01, 9 = 0.2, St = 1/128, At = 1, p = 0, 
q = 0.8. 
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